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The ground-state energies of 2, 3, 4 and 5 7r + 's in a spatial volume V ~ (2.5 fm) 3 are computed 
with lattice QCD. By eliminating the leading contribution from three-7r + interactions, particular 
combinations of these n-n + ground-state energies provide precise extractions of the tv + tt + scattering 
length in agreement with that obtained from calculations involving only two 7r + 's. The three- ir + 
interaction can be isolated by forming other combinations of the n-ix + ground-state energies. We 
find a result that is consistent with a repulsive three-7r + interaction for m w < 352 MeV. 



A major goal of strong-interaction physics is to deter- 
mine the spectrum and interactions of hadrons and nuclei 
from Quantum Chromodynamics (QCD). Lattice QCD 
is the only known way to rigorously compute strong- 
interaction quantities, and an increasing effort is being 
put into understanding the lattice QCD calculations that 
will be required to extract even the most basic properties 
of light nuclei. It is clear that at some level, the interac- 
tions among three or more hadrons play a significant role 
in nuclei, and an important goal for lattice practitioners 
is to determine the parameters of these interactions. We 
report on the first lattice QCD calculation of systems 
comprised of more than two hadrons. 

The simplest multi-hadron systems (both conceptu- 
ally, and from a numerical perspective) consist of n pseu- 
doscalar mesons of maximal isospin. Interactions among 
multiple pions are important to explore for phenomeno- 
logical reasons. Two- and three-pion interferometry is 
currently being used to determine the coherence of the 
pion source in heavy-ion collisions [T]. Further, multi- 
pion interactions impact the formation of a pion conden- 
sate which is energetically favored in systems with large 
isospin chemical potential, and will influence the proper- 
ties of (hot) pion gases. In this work we perform lattice 
QCD calculations of the ground-state energies of ir + ir + , 
7r + 7r + 7r + , 7r + 7r + Tr + 7r + , and 7r + 7r + 7r + 7r + 7r + systems in a 
spatial volume of V ~ (2.5 fm) 3 with periodic bound- 
ary conditions and a lattice spacing of b ~ 0.125 fm. 
These systems provide an ideal laboratory for investi- 
gating multi-particle interactions, as chiral symmetry 
guarantees relatively- weak interactions among pions, and 
multiplc-pion correlation functions computed with lattice 
QCD do not suffer from signal-to-noise issues that are ex- 
pected to plague analogous calculations in multi-baryon 
systems. The ir + ir + scattering length is extracted from 



the n > 2 pion systems with precision that is compa- 
rable to (and in some cases better than) the n = 2 de- 
termination (2]. Additionally, a result that is consistent 
with a repulsive three-pion interaction of magnitude ex- 
pected from naive dimensional analysis (NDA) is found 
for 7tv< 352 MeV. 

At finite volume, the ground state energy of a system 
of n bosons of mass M is shifted from its infinite vol- 
ume value, n M . In a periodic cubic spatial volume of 
periodicity L, this shift is known tobe[31ll[Sl[El[71IEllS] 
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where a and r are the two-boson scattering length and 
the effective range parameter, respectively, and 7/3 is the 
renormalization-group invariant (RGI) three-boson inter- 
action (773 is renormalization-scheme and scale indepen- 
dent, but depends logarithmically on L; in terms of the 
three-particle interaction defined in Ref [8J, rjjf = 773 (/w) + 
^^(3\/3-4 7 r)log(^ J L) - ^[2Q + TZ]). The geomet- 
ric constants appearing in Eq. ([IJ are X — —8.9136329, 
J = 16.532316 and K = 8.4019240. At this order, the 
energy is only sensitive to a combination of the effective 
range and scattering length, a = a + j^a 3 r and in what 
follows we replace a — > a, eliminating r. The above ex- 
pansion is valid provided a, r <C L with an additional 
constraint on n [Hj. 

Various combinations of the energy differences defined 
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in Eq. ([TJ) are particularly useful in what follows. First 
L 3 M(AE n m(m 2 - 3m + 2) - AE m n(n 2 - 3n + 2)) 
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x [-X 3 + 3XJ+(19 + 5mn-10(n + m))/C] j , (2) 

(for n, m > 2) is independent of 773 and allows a deter- 
mination of a up to 0(1/ L 4 ) corrections (combinations 
achieving the same result using all of the n = 3,4, 5 en- 
ergies can also be constructed). Second, the three body 
parameter can be directly determined from 
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(n > 2) with corrections arising at (D(l/L). Additionally, 
the dimensionless quantity 
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provides a useful check of the convergence of the expan- 
sion. Combinations involving A£? 2 ,3,5 and A£? 2j 4,5 that 
vanish at this order can also be constructed. 

The requisite ground-state energies are extracted from 
the n-TT + correlation functions defined by 



C n (t) 
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where x 7r +(x) = u a (x)jsd a (x) is an interpolating opera- 
tor for the 7r + (a is a color index). The sums in Eq. |5| 
project the correlation functions onto the represen- 
tation of the cubic symmetry group (in the continuum 
this corresponds to angular momentum £ — 0,4, . . .). As 
n increases, the number of Wick contractions involved 
in computing C n (t) increases as n\ 2 . In the limit of 
isospin symmetry, the correlation functions in Eq. |5| 
with n < 13 require the computation of only a single 
quark propagator, S(x; 0) (for n > 12 additional prop- 
agators are required to circumvent the Pauli exclusion). 
As an example, the n = 3 correlator can be expressed as 

C 3 (t) = tr [H] 3 - 3 tr [n] tr [II 2 ] + 2 tr [n 3 ] , (6) 

where II = ^ x 7sS'(x, t; 0)756'''' (x, t; 0) and the trace is 
over Dirac and color indices. 

In this work we have computed Ci 2,3,4,5 (£) in 
mixed-action lattice QCD, using domain- wall valence 
quark propagators from Gaussian smeared-sources on 



the rooted-staggered coarse MILC gauge-configurations 
(20 3 x 64) after HYP-smearing and chopping (see Refs. [2, 
ITU] for details). These are computed at pion masses of 
- 291, 352, 491, 591 MeV. Details of the propaga- 
tors used in the correlation functions are given in Table [TJ 
and can be found in Ref. [2]. 



m w (MeV) 


A'cfg 


iVsrc 




291.3(1.0)(1.0) 


468 


16 


1.990(11)(14) 


351.9(0.5) (0.2) 


769 


20 


2.3230(57)(30) 


491.4(0.4)(0.3) 


486 


24 


3.0585(49)(95) 


590.5(0.8)(0.2) 


564 


8 


3.4758(98)(60) 



TABLE I: Parameters of the domain-wall propagators used 
herein. A lattice spacing of 6 = 0.125 fm has been used to 
convert from lattice to physical units. The number of gauge 
configurations is N c i s , and the number of sources per config- 
uration is iV src . For further details see Ref. [2]. 



The energies of n pion states are dominated by the n 
single-pion energies, with the interactions contributing a 
small fraction of the total energy. To extract the resulting 
energy shifts, AE n , the ratios of correlators 



G n {t) 



C n (t) 



Ae 



-AE„ t 



(7) 



are formed, where the second relation holds in the limit of 
infinite temporal extent and infinite number of gauge con- 
figurations. Inclusion of the effects of temporal bound- 
aries (here Dirichlet boundary conditions are used) is 
complicated for multi-hadron systems, and our analysis is 
restricted to regions where an effective-mass plot clearly 
shows that the ground state is dominant. 

For the quantities discussed below, both jackknife and 
bootstrap analyses of the correlators and effective masses 
{e.g., log [G n (t) /G n (t — 1)]) are performed for each en- 
ergy or combination thereof. These are then used in 
correlated and uncorrelated fits to the t dependence to 
extract the relevant quantity. Our systematic uncertain- 
ties are determined by comparison of our different anal- 
ysis procedures and variation of the fit ranges. To avoid 
uncertainties arising from scale setting, we focus on the 
dimensionless quantities rn 7r a Tr + 7r + and rn^f^ff^ is 
expected to scale as m" 1 /" 4 by NDA). 

The 7r + 7r + scattering length (more precisely, the com- 
bination a 7r + 7r +) has been studied repeatedly in lattice 
QCD using the finite- volume formalism of Liischer [6] (for 
a/L <C la perturbative expansion gives the n = 2 case of 
Eq. 0). In particular, a precise extraction of this scat- 
tering length has been performed using the same prop- 
agators as used in this work [2j. Therefore, the utility 
of multi-pion energies in extracting the scattering length 
can be ascertained, as summarized in Figs. [TJ [2] and [3] 

In Fig. [TJ the energy shifts for n = 2, 3, 4, and 5 
are displayed for the highest-precision calculation, ~ 
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FIG. 1: Energy shifts, AE n , for n = 2, 3, 4, 5 tt+'s for the cal- 
culations with m-n ~ 352 MeV. The statistical and systematic 
uncertainties of the fits have been combined in quadrature. 



352 McV. Clear plateaus arc visible for each n; indeed, 
the relative uncertainty decreases with increasing n in the 
range explored (this is particularly clear for the calcula- 
tion with ~ 291 MeV). Since multiple combinations 
of pions interact in an n-pion state over a larger volume 
of the lattice, a statistically-improved signal results. 

Figure [2] presents extractions of the scattering length 
at all four orders in the 1 / L expansion in Eq. (JT|) for 
Tfiff ~ 352 MeV. For n > 2, the N 3 LO (l/£ 6 ) extraction 
is performed using Eq. ([2| with the point at n = 3 arising 
from the energy shifts AE4 and AE$, and so on. Sig- 
nificant dependence on n is observed in the lower-order 
extractions (LO, NLO and NNLO), indicating the pres- 
ence of residual finite- volume effects. However the most 
accurate extractions using Eq. which eliminates the 
three-7r + interaction (Eq. (JT|) for n = 2), are in close 
agreement for all n. This provides a non-trivial check of 
the n-dependence of Eq. (|T|), particularly the presence of 
a term that scales as | n I , which can be identified as the 
thrce-pion interaction. 
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FIG. 2: Extracted values of m 7r a 7r + 7r + at ~ 352 MeV. 
LO, NLO and NNLO correspond to extractions of a at 
0(l/L 3 ,l/L 4 ,l/L 5 ) from Eq. fi), respectively. The N 3 LO 
results for {n, m} = {3, 4}, {3, 5}7 and {4,5} are determined 
from Eq. |2|. For n — 2, the exact solution of the eigenvalue 
equation [of is denoted by "Liischer". 

The effective ra 7r a 7r + 7r + plots for {n,m} = {3,5} are 
shown in Fig. [3j Agreement is found at the level of cor- 



relation functions with those of n = 2, {n,m} — {3,4} 
and {n,m} = {4,5}. This agreement suggests that 
higher-order effects in X/L (such as higher-derivative in- 
teractions and four-particle interactions, which occur at 
0(1/ L 8 ) and 0(1/L 9 ), respectively) are small. For the 
calculations with ~ 291 MeV, the n > 3 effective 
plots are significantly "cleaner" than for n — 2. 
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FIG. 3: Effective m n a v + w + plot for {n,m} = {3,5} using 
Eq. pi). The statistical and systematic uncertainties of the 
fits have been combined in quadrature. 



To isolate the three-body interaction, we turn now to 
the combinations defined in Eq. (3]), and the effective 
m-nf^rjs plots are shown in Fig. [4 A nonzero value of 
m ^ftvz is found for ~ 291 and 352 MeV. Figure [5] 
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FIG. 4: Effective plots extracted from the n = 3, 4, 

and 5 7r + energy shifts. The fits shown correspond to the n = 
5 calculation. The statistical and systematic uncertainties 
have been combined in quadrature. 

and Table Ull summarize the results for the RGI three- ir + 
interaction, m^f^fj^, at L = 2.5 fm. The magnitude of 
the result is consistent with NDA. In Table [TTJ we also 
present f£ri3 (a* = V&)i a quantity that has a well- 
defined infinite- volume limit (unlike 77^) but is scale and 
scheme dependent. Its scale dependence is given below 
Eq. , and we use the MS subtraction scheme [5] . 
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FIG. 5: Mass dependence of the RGI three-7r + interaction 
m irftv3- The statistical and systematic uncertainties have 
been combined in quadrature. The vertical dashed line de- 
notes the physical value of m 7r + //„.+ . 



m n (MeV) 


291 


352 


491 


591 


W»/*7? 3 (/i = 1/6) 


1.3(2)(7) 
1.2(2)(7) 


0.8(1)(2) 
0.7(1)(2) 


0.4(2)(4) 
-0.1(2)(4) 


-0.4(3)(4) 
-1.3(3)(4) 



TABLE II: The 7r + 7r + 7r + interaction as defined in Eq. 
The most precise result (using n = 5) is quoted. 



Finally, Eq. Q and its counterparts involving other 
combinations of energies allow for a determination of 
residual 1 /L 7 contributions to the quantities we have ex- 
tracted at N 3 LO. They are all consistent with zero and 
are < 0.05. 

Further lattice QCD calculations are required before a 
definitive statement about the physical value of the three- 
pion interaction, m^f^fj^, can be made. While at lighter 
pion masses, there is evidence for a contribution to the 
various n-pion energies beyond two body scattering that 
scales as the three-body contribution in Eq. (JT|), a num- 
ber of systematic effects must be further investigated. 
The extraction of this quantity has corrections that are 
formally suppressed by a/L, however, the coefficient of 
the higher order term(s) may be large, and the next or- 
der term in the volume expansion needs to be computed 
(for n — 3, this result is known j9]). It is also possible 
that the signals seen in Fig. [4] are artifacts of the lattice 
discretization, but the observed scaling that is consistent 
with f n j suggests this is not the case. However, calcula- 
tions at a finer lattice spacings and with different lattice 
discretizations are required to resolve this issue. 

As the lattice QCD study of nuclei is an underlying 
motivation for this work, it is worth considering diffi- 
culties that will be encountered in generalizing the re- 
sult described here to baryonic systems. Certain difficul- 
ties have been discussed in Ref. [5]. Here we focus on 
the numerical issues. The ratio of signal-to-noise scales 
very poorly for baryonic observables |llj . requiring an 
exponentially-large number of configurations to extract 
a precise result. Also, the factorial growth of the com- 
binatoric factors involved in forming the correlators for 



large systems of bosons and fermions and the high powers 
to which propagators are raised {e.g., for the 12-7r + corre- 
lator, there is a term 43545600 tr[n 11 ]tr[n]) implies that 
the propagators used to form the correlation functions 
must be known to increasingly high precision. There is 
much room for theoretical advances in this area. 

In this work we have numerically studied the ground- 
state energies of n = 2, 3, 4, 5 7r + 's in a cubic volume 
with periodic boundary conditions using lattice QCD. We 
find that the ir + ir + scattering length can be extracted 
from combinations of these energies that eliminate the 
three-7r + interaction, and agree with previous n = 2 cal- 
culations [2 j . In some cases the precision of the extrac- 
tion is improved. We have found evidence of a repulsive 
three-7r + interaction for 352 MeV. Future calcula- 

tions will extend these results to larger n and to systems 
involving multiple kaons and pions. Further, calculations 
must be performed in different spatial volumes to deter- 
mine the leading correction (0(1/ L)) to the three- 7t + 
interaction, and at different lattice spacings in order to 
eliminate finite-lattice spacing effects, which are expected 
to be small. 
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